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The authors of a recent Letter [1] based their study of 
rogue waves in nonuniform currents on a modified non- 
linear Schrodinger equation (NLSE; see Eq.(l) in [1]). 
However, I will show below that equation is not correct. 
It gives wrong solutions even in the first order on the 
supposedly small parameter U/c s , where U(x) is a cur- 
rent, and c g = g/(2u>) [here u> is a mean frequency of 
a quasi-monochromatic wave train, and g is the gravity 
acceleration] . I also suggest an accurate variant of NLSE 
[see Eq.© below], valid in the presence of a large-scale 
nonuniform current under condition (1 + AujU/g) > 0.2. 

The key point in the derivation is that with non- 
small U/cg we may not assume a globally narrow wave- 
number spectrum, since a local wave-number k essen- 
tially [up to several times!] varies along x together with 
U [see Eq.([2|) below]. First, we should find monochro- 
matic solutions of the following linearized water-wave 
problem: r) t = -{Urj) x + kip, and -ip t = Uip x + grj, 
where rj(x,t) is a vertical displacement of the free wa- 
ter surface from a steady-state profile, ip(x,t) is a sur- 
face value of the wave velocity potential, and a linear 
pseudo-differential operator k multiplies the Fourier im- 
age ijjk{t) by \k\. Assuming i] — Re[Q(a;, oj) exp(— iiot)} 
and ip = Ke[P(x, uS) exp(— iu;t)], we have 

u?P + iuj(d x U + Ud x )P - d x U 2 d x P = gkP. (1) 

Since U (x) changes over many typical wave-lengths, we 
can write an approximate solution of Eq. ([l} in the form 

P(x,u) « ^(x, uj) exp if x k(ui,U)dx , where ^(x,ui) 
is a slowly varying function, and a positive function 
k(uj,U) satisfies the local dispersion relation for deep- 
water waves, u> = Uk + \fgk. In an explicit form, 

k(ui, U) = [g + 2u)U - y/g 2 + Agu)U]/{2U 2 ). (2) 

It is important that with positive k we have kP w —id x P. 
Substitution into Eq.([T]) gives us the equation 

uj{U x ^ + 2U^ x )-2kU(y x ^ + U^ x )-U 2 k x ^ + g^ x w 0. 

Multiplying it by ^ and integrating on x, we obtain [ojU— 
kU 2 + g / 2]fy 2 ps const (it is equivalent to the wave action 
conservation [2]), or <3> « -iC[l + 4ujU/g]~ 1/4 , with a 
complex constant C. Since gQ = iLuP — UP X , we have 
Q(x, u) = CM(x, u) exp [i J x k(u, U)dx\ , where M w 
(fc/g) 1 ^ 2 [l + iujU/g]^ 1 ^ 4 . What is important, at small 
U the wave amplitude behaves as M/Mq ps (1 — 2ujU /g), 


while in [1] the corresponding factor is exp[— U/(2c s )] ps 
(1 — coU/g), which is not correct. 

Let us now consider a linear superposition of the mono- 
chromatic solutions in a narrow frequency range near w, 

v = RcJ d£C(£)M(x, + e {-^+t)t+* F Hu,+t,u)dx} 
w Re{Q[x,t)M(x,u)e{- iut+i S Xk{ - u W dx }], (3) 
where Q(x, t) is defined by the following integral, 


Since the frequency spectrum C(£) is concentrated at 
small £ in a range AQ <C w, we can expand [k(uj + £, U) — 
k(u>, U)] in powers of £, and thus derive a partial-differ- 
ential equation for <d(x,t) in the linear regime, — iO x = 
ik u Qt — 0-/2)k uu &tt + • • •■ Here k^ and fc ww are partial 
derivatives of k (w, U): k u = [4w/g)[l + v + (1 + v) 1 / 2 ]" 1 
and kuu, — (2/g)(l + u) -3 / 2 , with v = 4uU(x)/g [wave 
blocking takes place at u* = — 1]. 

Eq.([3]) means that the wave envelope is A(x,t) ps 
Q(x,t)M(x, u>) [A here is not the same as A in [1], but 
\A\ = \A\]. It results in the following estimate for the 
quantity I = (k\A\)(k/ AK), where AK m k^AQ is a 
width of the local spectral fc-distribution, 

I v /I = 2%[l+v/2+Vl + v]- r t{Vl + v+l)(l+v)^. (4) 

Thus, Eq.(8) in [1] should be corrected as written below, 


A max (v)M /(M v A ) 


l + 2Jl-[V2eNI v /I }- 2 . (5) 


However, applicability of formula ((H) implies (Ail) 2 <C 
6fc w / 'kuw', otherwise AK s6 k^Afl. Since practically im- 
portant are values Aft ps (0.1...0.2)w, Eq.Q can be good 
at v > —0.8 only. The same condition arises when we 
require that the neglected linear higher-order dispersive 
terms are small comparatively to k U u®tt- 

To complete our derivation of NLSE with variable co- 
efficients for a weakly nonlinear deep-water wave train 
in a large-scale nonuniform current, we have in a stan- 
dard manner to take into account the nonlinear fre- 
quency shift, which for fixed k is well known to be 
Sui ps \fgkk 2 \A\ 2 /2. For fixed ui, that corresponds to 
the nonlinear wave-number shift 8k ss —k u yfgkk 2 \A\ 2 /2. 
Using the relation \A\ ss |0|M, we finally derive 

1 k k 3 \/k 
t Q x+l k^Q t --k^Q u --^=y=\Q\ 2 e^0. (6) 

2 2y/g + AojU 
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